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I.   INTRODUCTION 

The  classic  single  period  inventory  problem  is  called  by 
various  well-known  names  such  as  "Newsboy  problem,"  "News- 
paperboy  problem,"  and  so  on.  [6].   In  this  thesis  it  will 
be  called  the  Newsboy  problem. 

Solutions  to  the  Newsboy  problem  are  well  known  under 
both  "risk"  and  "uncertainty"  conditions.  [7].   The  purpose 
of  this  thesis  is  to  study  the  Newsboy  problem  under  initial 
conditions  of  "Uncertainty,"  (i.e.,  when  the  possible  future 
for  demand  is  not  yet  known)  and  then,  for  the  subsequent 
periods  as  data  gathers,  to  suggest,  structure,  and  evaluate 
decision  procedures  during  this  transition  phase  to  decision 
under  "Risk." 

The  inventory  model  studied  in  this  thesis  is  a  time 
independent  (stationary)  and  linear-cost  type.   The  study 
will  consider  the  problem  of  minimizing  the  expected  cost. 
Since  with  the  Newsboy  problem  both  the  expected  profit 
maximization  and  expected  cost  minimization  cases  have 
identical  solutions,  the  study  should  be  applicable  to  both 
interpretations  of  the  problem. 

Chapter  II  reviews  the  Newsboy  problem  for  decisions 
under  risk  and  decisions  under  uncertainty,  and  then  suggests 
possible  alternative  procedures  for  "from  uncertainty  to 
risk"  conditions.   Among  these,  Order  Statistics  (a  non- 
parametric  approach)  appear  promising  and  are  discussed  in 


detail  in  Chapter  III.   Chapter  IV  evaluates  the  alternative 
nonparametric  transition  phase  procedures  of  the  Newsboy- 
problem  and  as  a  result  proposes  a  transition  phase  pro- 
cedure.  In  Chapter  V  conclusions  and  recommendations  for 
further  work  are  given. 


II.   SOLUTION  OF  THE  SINGLE -PERIOD  INVENTORY  (NEWSBOY) 
PROBLEM  UNDER  RISK  AND  UNCERTAINTY  CONDITIONS 


In  this  section  the  Newsboy  problem  is  examined  under 
risk  and  uncertainty,  and  optimal  decision  rules  for  these 
conditions  are  given.   These  two  cases  form  the  bounds  of 
the  transition  phase  and  permit  us  to  characterize  the 
transition  phase  decision  problem  we  intend  to  study. 

Here,  when  it  is  said,  "under  risk,"  it  is  meant  that 
the  probability  distribution  of  demand  is  known.   When  the 
distribution  of  demand  is  not  known,  and  data  in   hand  is 
not  sufficient  to  estimate  the  possible  futures  of  demand, 
then  the  problem  is  called  "under  uncertainty."  [7], 

A.   SOLUTION  OF  THE  NEWSBOY  PROBLEM  UNDER  RISK 

Let  fr>(x)  be  the  probability  density  function  of  the 
demand  for  continuous  populations,  and  pn(x)  be  the  proba- 
bility mass  function  for  discrete  populations.   If  C    is 
the  stock-out  or  shortage  cost  per  unit  of  unsatisfied 
demand,  and  C   the  carrying  or  surplus  cost  per  unit  of 
excess  supply,  then  the  cost  for  a  particular  period  of 
Newsboy  problem  will  be: 

(CC(S-X) ,  if  0  <  X  <  S,  and 

cost  ccs)  =jcso(x-S),  if  x  >  S, 

where  X  is  the  realization  of  demand  for  that  period,  and 
S  is  the  inventory  level  in  hand.   Here  S  represents  the 
decision  variable  for  this  problem. 
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The  expected  cost   E[C(S)]   which  results  from  keeping 
an  inventory  level  S  in  hand,  will  be: 


S 
E[C(S)]  =  Z    Cc(S-X)pn(X)  +  S   C   fX-S)   p  (X)        (1) 
X=0        U  X=S+lso        u 


for  discrete  demand  distributions  and, 

S 
E[C(S)]  =  Cc  /(S-x)  fD(x)dx  +  CsQ  /  (x-S)   fD(x)dx,   (2) 

0  S 

for  continuous  demand  distributions. 

Now,  the  expected  cost  for  a  period  is  to  be  minimized 
and  we  wish  to  find  an  optimal  value  of  S,  say  S  ,  that  will 
yield  the  minimum  expected  cost.   In  other  words,  (1)  and 
(2)  are  to  be  minimized  with  respect  to  S. 

In  the  discrete  demand  distribution  case,  two  necessary 
conditions  for  a  minimum  at  S   are: 

E[C(So+  1)]  -  E[C(SQ)]  >  0,  and 

B[C(S0-1)]  -  E[C(S0)]   >  0. 

These  conditions  yield  the  well  known  result  [7]: 

as  conditions  for  a  minimum  expected-cost  solution  at  S  . 

r  o 

Similarly,  for  continuous  demand  distributions,  to 
minimize  (2),  the  first  derivative  of  E[C(S)]  is  taken  with 
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respect  to  S,  and  after  it  is  set  equal  to  zero,  yields 

C 


so 


c     SO 


where  Fn(S  )  is  the  cumulative  distribution  function  (CDF) 
of  demand  at  S  .   Thus,  the  S  values  which  satisfy  equations 

0  0  /    -i 

(3)  and  (4)  are  the  optimum  inventory  levels  to  keep  on  hand 
in  order  to  minimize  the  expected  cost. 

This  solution  (under  risk)  is  quite  idealistic,  since  it 
assumes  that  the  distribution  of  demand  is  known.   In  fact, 
in  reality  we  are  seldom  certain  about  the  distribution.   In 
most  cases,  neither  the  distribution  of  demand,  nor  an  ade- 
quate amount  of  data  to  provide  a  good  basis  for  estimation 
of  the  distribution  is  initially  available.   The  situation 
when: 

1.  Distribution  of  demand  is  not  known,  and 

2.  Enough  data  from  the  system  to  estimate  the  distri- 
bution of  demand  is  not  available,  but 

3.  We  do  have  some  idea  about  the  range  of  demand  D, 

say  0  <  D  <  D    ,  is  one  we  shall  call  "uncertainty." 

J  max  ' 

B.   SOLUTION  OF  NEWSBOY  PROBLEM  UNDER  UNCERTAINTY 

It  has  been  shown  that  whether  minimax  cost,  minimax 
regret,  or  the  Laplace  approach,  assuming  a  uniform  distri- 
bution of  demand  is  chosen  as  a  principle  of  choice,  the 
optimal  decision  under  uncertainty  is  the  same  [6],  [7], 
namely, 
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C     D 
_       so   max  r  j.  •      j  „   j 

S   =   p p ,  for  continuous  demand, 

0     Lc  so 

S   <    Cso(Dmax  +  ^     s    , 
c    so 


for  a  discrete  demand  distribution.   The  results  of  these 
uncertainty  procedures  are  not  affected  by  newly  generated 
demand  data  unless  it  causes  one  to  revise  the  range  of  demand. 

C.   TRANSITION  PHASE  OF  THE  NEWSBOY  PROBLEM 

Now  the  two  bounds  or  limit  solutions  (uncertainty  and 
risk  conditions)  for  the  transition  phase  of  the  Newsboy 
problem  are  reviewed.   Implicitly,  to  proceed  in  accordance 
with  these  current  solutions,  one  begins  under  Uncertainty 
and  uses  the  Uncertainty  optimum  inventory  level  S   (based 
on  the  range  of  demand) ,  until  demand  data  is  adequate  or 
complete  to  estimate  the  distribution  of  demand  and  change 

to  the  optimum  inventory  level  S   for  the  Risk  conditions. 

r  J  o 

As  mentioned  before,  the  purpose  of  this  study  is  to  find  a 
good  transition  phase  procedure  for  the  intermediate  periods 
when  some  data  has  accumulated  and  is  available.   A  transition 
phase  procedure  may  be  characterized  by  two  decisions: 

1.  When  to  begin  the  transition  phase  (i.e.,  at  what 

period  should  we  leave  Uncertainty  procedures  to  begin  to 

use  an  estimator  of  S  ) . 

o 

2.  What  estimator  of  S   to  use. 

o 

After  explaining  the  nature  of  the  transition  phase  in 
the  next  chapter,  three  candidate  estimators  for  the  optimum 
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inventory  level  S   to  use  during  transition  phase  will  be 
proposed. 

As  we  remember,  the  optimal  choice  of  S  under  Risk  was 
based  on  a  quantile.   Let  us  call, 


Cso 

P  =  C   +  C 

c    so 


and  then  from  equation  (4)  we  have, 


S   =  X   =  F"1  (p)  , 
o    p     D   ^ 


which  provides  a  minimum  expected-cost  solution  under  Risk. 
Thus,  the  problem  is  to  locate  the  p   quantile  of  the 
Cumulative  Distribution  Function,  Fn(X)  .   As  more  and  more 
data  is  collected,  one  might  be  able  to  estimate, 

-».  /\  /N  _  "I 

S   =  X   =  F   (p)  , 
o    p       r 

as  a  function  of  the  data.   Some  alternative  approaches  to 
estimate  this  quantile  were  considered  to  be: 

1.  To  consider  similar  systems  for  which  we  have  demand 
data,  and  then  to  apply  a  Bayesian  approach. 

2.  To  decide  about  the  form  of  the  distribution  of 
demand  and  carry  out  a  maximum  likelihood  procedure  to  find 
the  parameters. 

3.  To  apply  a  nonparametric  approach. 

For  the  first  two  cases,  if  one  is  able  to  decide  about 
the  form  of  distribution  such  as  lognormal,  gamma,  compound - 
poisson  process,  beta,  etc.,  then  the  p   quantile  of  demand 
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might  be  estimated  quite  well  subject  to  data.   It  is  the 
third  approach,  that  of  nonparametric  statistics,  which  is 
the  center  of  interest  in  this  thesis,  and  forms  the  basis 
for  the  chapter  that  follows. 
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III.   A  NONPARAMETRIC  APPROACH  FOR  THE  NEWSBOY 
PROBLEM,  TRANSITION  PHASE 

In  this  section  the  possibilities  for  applying  a  non- 
parametric  or  distribution-free  approach  to  the  problem  of 
the  transition  phase  of  the  Newsboy  problem  will  be  explored. 
After  explaining  the  reasoning  for  a  nonparametric  approach, 
three  candidate  estimators  of  the  optimum  inventory  level 

S   for  X  )  will  be  proposed, 
op 

A.   WHY  A  NONPARAMETRIC  APPROACH? 

As  mentioned  at  the  end  of  the  preceding  section,  if  the 
decision  maker  is  able  to  predict  the  form  of  the  demand 
distribution,  then  a  parametric  approach  can  be  more  efficient 
and  have  a  higher  convergence  rate  to  the  true  value  of  the 
estimator  than  a  nonparametric  equivalent.   Hadley  and 
Whitin  [5] ,  give  various  considerations  relating  to  choice 
of  a  demand  distribution  and  about  inventory  demand  predic- 
tions, in  detail.   For  a  parametric  approach,  such  as  would 
result  from  fitting  a  theoretical  distribution  to  data,  they 
suggest  to  use  poisson  distribution  when  the  lead  time  is 
relatively  low,  or  normal  distribution  when,  it  is  not  rela- 
tively low,  etc. 

Since  it  is  quite  difficult  to  specify  the  form  of  the 
demand  distribution,  a  nonparametric  approach  was  considered 
to  be  of  particular  interest  for  this  problem.   Also,  Hadley 
and  Whitin  point  out  that,  for  single-period  inventory  models, 
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like  the  Newsboy  problem,  the  empirical  distribution  function 
(nonparametric)  is  reasonably  applicable,  unless  we  are  forced 
to  operate  the  model  at  the  tail  of  the  distribution.  [5]. 
Performance  of  the  model  decision  rules  was  assumed  to 
be  critical  during,  say,  the  first  50  periods  when  data  is 
inadequate.   The  decision  maker  needs  to  operate  the  system 
well  during  that  critical  time.   Later  on  he  will  have  a 
sufficient  amount  of  data  to  treat  the  problem  as  one  under 
Risk  and  simply  use  the  optimal  decision  rules  for  Risk.   For 
these  early  periods  (transition  phase)  a  nonparametric  approach 
was  used. 

B.   THE  EFFECT  OF  THE  ESTIMATORS  OF  OPTIMUM 
INVENTORY  LEVEL  ON  EXPECTED  COST 

Before  we  suggest  some  estimators  of  optimal  inventory 
level,  X  ,  it  may  be  useful  to  see  how  those  estimators 
affect  the  total  expected  cost  in  a  period.   In  the  work 
that  follows,  it  is  assumed  that  the  demand  process  is  time 
independent  (stationary) ,  which  means  that  seasonal  effects 
and  trends  are  removed.   Also,  we  assume  that  for  the  previous 
period  we  know  what  the  demand  was. 

Before  we  continue  on  with  the  mathematical  development, 
it  is  useful  to  expand  on  the  conditions  for  an  optimal 
expected-cost  solution.   When  the  demand  distribution  is 
assumed  known,  we  remember  from  equation  (2)  that  the  expected 
cost  for  the  case  of  continuous  demand  is: 


E[C(S)]  -  Cc  /  (S-x)   fD(x)dx  +  CsQ  /  (x  -  S)   fD(x)dx. 
o  S 
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To  find  an  optimum  inventory  level  S  ,  taking 


b-[ccs)]  -Aimm 


dS 


S  =  S 


(5) 


and  setting  it  equal  to  zero  (necessary  conditions  for 
minimization) ,  we  obtain: 


C    +  C     Fn(SJ  +  C„    Fn(SJ  =  0  ,  or 
so    so    Dv  oJ  c    D  o^      ' 


so 


FDtSo)  ~  C_  -  C 


so 


=  p  ,   and   SQ  =  FD  (p) 


Similarly, 


E~[C(S)]  =  d  ElC(S^ 


d  S' 


S  =  S 


•  c   fnrs  )  +  c   fn(s  ) 

so    D^    oJ         c    D^  oJ 


=  (C   +  C  )   fnCSJ  >  0  , 
v  so    c     D  o 


(6) 


which  forms  a  sufficient  condition  for  minimization.   Here, 

S   is  the  value  that  corresponds  to  the  p    quantile  of  the 

population,  say  X  . 

We  are  interested  in  the  case  where  the  population  is 

not  known,  but  there  exists  a  sample  of  size  n  of  demand 

data,  acquired  during  the  last  n  periods. 

Let  Xt  ,  X-,,....,  X  be  realization  of  demand  during  the 
l   z       n 

last  n  periods.   Although  we  don't  know  the  value  of  X  yet, 
we  do  have  a  chance  to  estimate  this  value  as  a  function  of 
available  data,  say: 


-1 


X  =  gCX1,X2,....,Xn)  =  FD"  (p)  = 


(7) 
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To  express  the  expected  cost  as  the  function  of  X  ,  let 


<5  =  S   -  S   =  X   -  X 
oo    p    p 


Equation  (2)  then  becomes: 


S  +  6 
E[C(SQ)]  =  Cc  f°    (SQ+6-x)   fD(x)dx  +  CS0 


+  C     /  (x  -  S  -  3)  fn(x)dx 

SO      g   +5       O  U 

o 


(8) 


If  equation  (8)  is  expanded  using  the  Taylor  series  around 
SQ,  it  gives 

E[C(SQ)]  =  E[C(So)]  +  E'[C(So)]   E[6]  +  E~[C(S0)]   ^-1  + 


Since  E'[C(S  )]  =  0  by  equation  (5),  and  having  E"[C(S  )] 
from  equation  (6) ,  then  for  three  terms  the  Taylor  expansion 
is : 


E[C(S  )]  =  E[C(S  )]  +  (C   +C  )   fn(S  )   Et ,  3 
L  "■  o  o      v  so    c     D   o     2! 


(9) 


Equation  (9)  shows  that: 

1.  Under  Uncertainty  conditions  the  expected  cost  for 
a  period  is  linear  in  the  mean  square  error  of  the  p 
quantile  estimators,  i.e.,  in 

E[52]  =  E[(Xp  -  Xp)2] 

2.  If  one  is  able  to  find  a  consistent  estimator. 

X  (n) ,  of  X  ,  then  when  n  gets  large  the  mean  square  error 
will  converge  to  zero   and  E[C(S  )]  will  converge  to  E[C(S  )], 
providing  that  the  moments  of  X  (n)  exist  and  behave  properly. 
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3.   The  optimal  performance  for  an  estimator  here  is  to 

yield  the  expected  total  cost  under  Risk,  E[C(S  )] . 

Now,  we  wish  to  find  a  good  consistent  point  estimator 

X   for  X  .   Also,  it  is  not  required,  but  it  is  desirable  to 
P      P 

obtain  one  that  is  an  unbiased  estimator  of  X  .   If  so, 

EKX   -  X  )  1  will  be  equal  to  the  variance  of  estimator  X  . 
L   p    P  P 

Then,  by  equation  (9),  if  it  is  possible  to  improve  this 

estimator  reducing  its  variance,  it  will  directly  reflect 

and  will  decrease  the  expected  cost  accrued  by  using  that 

estimator. 

To  obtain  candidate  nonparametric  estimators  of  X  , 

order  statistics  will  be  used.   In  the  following  section 

we  shall  develop  three  different  extimators  of  X  ,  using 

an  order  statistics  approach.   Additional  details  about  order 

statistics  and  empirical  distribution  functions  related  to 

quantile  estimation  are  given  in  Appendix  A. 

1 .   The  First  Candidate  Estimator  of 
Optimum  Inventory  Level:   X,  , 

Let  Xf1N  <  Xro.  <  .  .  .  <  X,  -.  ,  be  a  set  of  n  reali- 
Cl)     (2J  (nj 

zations  from  an  unknown  distribution  in  order  statistics 
form,  where  X,  ,  is  the  r   largest  value  (r   order 
statistic) .   Let  r  be  defined  in  the  following  way: 

n  p  ,  when  n  p  is  integer 
r  = 


[np   +1]  ,.  otherwise, 

where  p  is  the  quantile,  n  is  the  number  of  periods  so  far, 

and  [X]  denotes  the  integer  part  of  X.   For  large  n,  X,  , 

is  a  consistant  and  unbiased  point  estimator  of  X   (Appendix  A) 

19 


In  Figure  1,  the  r  order  statistic  X,  •>  is  represented 
graphically  with  the  distribution  function  F„(X)  and  the 
empirical  distribution  function  G„(X) . 

Know  results  of  this  estimator  include  the  following. 
For  a  uniform  (0,1)  distribution, 

E[X(r)]  =  E[U(r)]  =  HTT  '  (10) 

and 

cov[Xr    VX        ]    =   cov[Uf   vUf    J    =   r(n-s+l) ,    for   r   <    s 

m  ls)  Lrj   lsJ     (n+l)^(n+2) 

(11) 

Similarly,  for  any  other  distribution, 

EtX(r)]  *  FD_1^  ,  (12) 

and 

cov[x      x     n  r(n-52+1)      •        .1\ n-^— 

™       ts)    (n+l)2(n+2)    fD(FD1(KIT)).fD(FD1(ir|T)) 

provided   r  <  s.  *■  ' 

Later  these  results  will  be  used  to  compare  the 

efficiency  of  the  estimators  of  X  . 
1  P 

2.   The  Second  Candidate  Estimator 

of  Optimum  Inventory  Level:   X,  >. 

As  a  second  candidate  estimator  of  X  ,  X,  .  is 

P    (r) 

offered  to  smooth  the  empirical  distribution  function  by  a 
linear  interpolation.   It  might  be  considered  that  the  area 
between  Gn(X)  and  the  actual  distribution  function  F„(X) 
may  be  a  measure  of  deviation.   So,  here  interpolation  was 
done  at  the  midrange,  yielding  a  new  empirical  distribution 
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function  R,  which  would  possibly  have  a  smaller  area  between 

R  and  FD(X),  than  the  area  between  GD(X)  and  F  (X)  .   A  possible 

realization  for  X--  -.  is  shown  in  Figure  2.   Since  the  area 

between  those  distribution  functions  is  subject  to  change 

because  of  randomness,  the  results  for  that  case  would  be 

evaluated  numerically  rather  than  analytically. 

3 .   The  Third  Candidate  Estimator 

of  Optimum  Inventory  Level:   X,  -. 

As  a  third  candidate  estimator  of  X  ,  a  pooled  order 

P    F 

statistic,  combining  other  order  statistics,  is  proposed. 

For  such  an  approach  to  be  successful  by  equation  (9) ,  it 

has  to  reduce  the  mean  square  error  to  a  lower  value  than 

the  r   order  statistic  does,  if  it  is  possible  to  make  this 

pooled  estimator  an  unbiased  estimator  of  X  .   Then  the 

P 

problem  is  one  of  comparing  the  variance  of  this  estimator 
with  the  variance  of  the  r   order  statistic. 

To  examine  this  subject,  it  is  reasonable  to  begin 
with  the  uniform  (0,1)  distribution,  since  the  expectations, 
variances,  etc.  of  its  order  statistics  are  simple  and 
exact,  and  then  to  discuss  the  assumptions  that  are  necessary 
before  we  generalize  it  for  other  distributions. 

a.   Xf.  for  the  Uniform  Distribution 

Let  X,   ,.,  X,  -v,  Xf   i-)  be  the  order  statistics 
from  the  uniform  (0,1)  distribution.   In  addition  to  the  r 


order  statistic,  a  pooled  estimator  in  the  form  of  a  linear 
combination  of  the  (r-1),  (r) ,  (r+1)  order  statistics  will 
be  considered.   Here  cur  third  candidate  estimator  is: 
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X(r)  =  Xp  =  a   X(r-1)  +  6   X(r)  +  *   X(r  +  1)  ■ 
where   a  +  g  +  .y  *  1. 

As  mentioned  earlier,  if  this  estimator  can  be  an  unbiased 
estimator  of  X  ,  then  a  comparison  of  the  variances  of  the 
r   order  statistic  (First  candidate  estimator) ,  and  this 
pooled  estimator  (Third  candidate)  will  indicate  which  one 
has  better  performance. 

To  maintain  unbiasedness ,  we  must  have: 


E[X(r)]  =  E[X(r)]  ,    or 

E[f,r)]  =  o.      E[X(r_i:)]  +  3   E[X(r)]  +  Y   E[X(r+1)] 


If  demand  has  a  uniform  distribution,  then 

E[f(r)l  -  a     C&)  *   B   C^)  ♦  r  eg*)  . 

This  quantity  is  equal  to  E[X,  .]  =  —pr   only  when  a  =  y. 
One  way  to  decide  on  the  values  of  a,  8  and  y    is  to  choose 
values  so  that  the  variance  is  minimized.   So,  we  would  like 
to  minimize 


Var(X(r))  =  Var(a   X(r_1}  +  B  X(r)  +  y   X(r+1))  , 
subject  to, 

a  +  6  +  y  =  1       and    a  =  y 

By  the  constraints  on  a,  0  and  y,  these  coefficients  can  be 
written  in  terms  of  a  only,  yielding 


Var 


(X(r))  =  Var{a(X(r.1)  *  X(r+1))  *■  (l-2a)  X(r)} 
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or, 

VarCf(r))  =  a2CVar(XCr_i:))  +  Var (X(r+1) ) ) 

♦  2  CovCX(T.1)fx(r+1)D 

+  (2a-4a2)(Cov(X(r_i:),X(r))  +  CovfX^  ,X  (r  +  1) ) ) 

+  (1  -  4a  +  4a2)Var(X(-r))  .. 

Taking  the  derivative  of  Var  (X,-.)  with  respect  to  a  and  setting 
it  equal  to  zero  gives, 

a  =     2  Var(XCr])  -  CovCX^.^  ,X(r))  -  Cov(Xm  ,Xfr  +  n) 
VarCXCr.1})  +Var(X(r+1))  *  2Cov(X(r.13  ,X(r+1) )  - 

4[Cov(Xr.t),XCr))  +Cov(X(r),X(r  +  i:))  -  Var(X(r))] 

(14) 
Having  available  the  values  of  variance  and  covariances  by 

equation  (11) ,  equation  (14)  becomes 

2r(n-r+l)  -  (r-l)(n-r+l)  -  r(n-r) 

a  =  

(r-1) (n-r+2)  +  4r(n-r+l)  +  (r+1) (n-r)  -  4(r-l)  (n-r+1)  - 

4r£n-r)  +  2(r-l) (n-r) 

the  n,r  terms  in  the  numerator  and  in  the  denominator  cancel, 

resulting  in  a  =  1/2.   So,  to  minimize  the  variance  of  L  -.  , 

the  coefficients  have  to  be   a  =  y  =  1/2,  and  3=0. 

This  says  that  when  a  population  is  uniform,  a  third 

reasonable  estimator  of  X   is: 

P 


X 


(r)  =  ^V-D  +  X(r+1))  »  <15> 
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yielding 


Var(X 


(r) 


2r(n-r+l)  -  (n+1) 
2(n+l)2  (n+2) 


If  the  ratio  of  the  two  candidate  estimators  is  taken, 


Var(X(r))       2r(n-r+l) 
VarCl^)      2r(n-r+l)  -  (n+1) 


>   1  , 


always.   This  implies  that  for  a  uniform  distribution, 
Var  (X,  O  >  Var  (X(  •«),  and  Xf  «.  is  the  better  estimator  of 
X  .   In  Figure  3,  for  n  -  20,  these  variances  are  sketched. 
When  r  is  close  to  bounds  (either  close  to  1  or  n) ,  the 
difference  between  those  two  variances  become  larger. 


(Variance) 
0.015 


Var(X(r)) 


0.002  ■• 


(r) 


FIGURE  3.   Variances  of  X.  »  and  X,  »  for  Uniform(0,l) 
Distribution, wnen  n  is  twenty  . 
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These  results  apply,  of  course,  only  to  uniform 
demand  distribution.   In  the  following  paragraphs  the  gener- 
alizations to  any  demand  distribution  will  be  discussed. 

b.   Generalization  of  X,-  -.  to  any  Demand  Distribution 

For  distributions  other  than  uniform,  it  is 
difficult  to  claim  that  X ,  , ,  as  given  by  equation  (15),  is 
exactly  unbiased.   Also,  the  variance  and  covariances,  as 
equation  (13),  will  be  functions  of  the  density  f^fX).   Even 
though  these  two  problems  exist,  some  reasonable  approxima- 
tions can  be  discussed.   If  we  have  both  a  large  n  and  fairly 
flat  densities,  since  (r-1),  (r) ,  (r+1)  are  the  closest 
neighbors  among  the  other  order  statistic  values,  the  densi- 
ties corresponding  to  these  points  would  tend  to  equal  each 
other  and  during  the  calculation  of  a,  as  in  equation  (14) , 
would  essentially  cancel. 

If  we  have  flat  densities  (where  the  three  points 
are  nearly  linear) ,  or  if  we  have  a  large  number  of  observa- 
tions (when  the  points  are  close  together  accordingly), 
X--  ->  would  be  approximately  unbiased,  too.   Thus,  it  may  be 
reasonable  in  some  cases  to  use  X<-  ->  for  other  distributions 
besides  the  uniform. 

Incidentally,  we  would  not  like  to  try  another 


pooled  estimator  of  X   such  as 

P 


X(r)  =  a2   X(r=2)+al  X(r-1)+6  X(r)+^   X(r+D  +  Y2  X(r+2) 
For  this  estimator,  the  assumptions  which  were  made  to  permit 
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X.  .  tcapply  the  distributions  other  than  the  uniform  may 
no -longer  be  appropriate.   For  example,  in  equation  (13), 

^^D1  Cn5T^  and   £D^FD   ^nTT^  would  Probably  have  quite 
different  values,  because  these  are  the  densities  of  the  two 
farthest  order  statistics,  X, r.2>  anc^  xrr+2")  *n  ^frV 
Otherwise,  if  these  densities  were  known,  we  would  not  use 
a  nonparametric  approach. 

So  far,  solutions  of  the  Newsboy  problem  under 
"Uncertainty"  and  under  "Risk"  have  been  reviewed,  and 
possible  procedures  for  a  transition  phase  between  them 
discussed.   Also  in  this  chapter,  three  candidate  estimators 
of  the  optimal  inventory  level  S   (or  X  )  were  developed  for 
use  as  alternative  transition  phase  procedures.   In  the  next 
chapter,  the  performance  of  these  alternative  transition 
procedures  will  be  evaluated  by  a  computer  simulation  of 
Newsboy  problem  decision  making.   Finally,  from  the  simula- 
tion we  will  obtain  information  about  what  time  to  leave  the 
minimax  (Uncertainty)  rule,  for  each  transition  phase  pro- 
cedure. 
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IV.   EVALUATION  OF  THE  CANDIDATE  ESTIMATORS 
OF  OPTIMUM  INVENTORY  LEVEL,  FOR  THE 
TRANSITION  PHASE  OF  THE  NEWSBOY  PROBLEM 


In  Chapter  III,  the  rationale  for  using  a  nonparametric 

approach  for  the  transition  phase  of  the  Newsboy  problem  was 

discussed.   It  was  shown  that  under  Uncertainty  the  expected 

cost  for  a  period  was  approximately  linear  in  the  mean  square 

error  of  the  optimum  inventory  level  S   (or  X  )  corresponding 

to  the  p   quantile,  and  three  candidate  estimators  of  X   or 

S  were  introduced, 
o 

In  this  chapter,  the  performance  of  these  estimators 
will  be  tested  for  several  demand  distributions  by  simulation, 
After  explaining  our  reasons  for  using  simulation,  we  will 
explain  the  structure  of  the  simulation  program  used.   In 
part  C,  the  simulation  parameters  used  in  the  program  are 
described,  together  with  the  reasons  for  choosing  these 
parameters.   Part  D  illustrates  the  results  of  simulation. 
These  results  will  be  used  later  to  compare  the  performance 
of  the  estimators  of  X  and  the  transition  phase  procedures, 
and  to  see  whether  there  were  any  apparent  adverse  effects 
of  approximations  or  assumptions. 

A.   RATIONALE  FOR  SIMULATION 

The  purpose  of  the  study  was  to  find  a  transition  phase 

procedure  (i.e.,  which  estimator  of  X  or  S   to  use  and  when) 

c  *■     '  p     o  ' 

during  the  transition  from  Uncertainty  to  conditions  under 
Risk  for  the  Newsboy  problem.  As  explained  before,  the 
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performance  of  an  alternative  transition  phase  procedure 
depends  directly  on  the  performance  of  the  candidate  esti- 
mator of  X  used,  in  that  there  is  a  one-to-one  correspondence, 
p 

A  primary  reason  for  simulation,  of  course,  was  to  study 
the  question  of  how  long  one  should  continue  using  the  mini- 
max  procedure  (for  Uncertainty)  before  changing  to  a  transi- 
tion phase  using  a  quantile  estimate.   Simulation  over  several 
demand  distributions  should  provide  insight  regarding  when 
(with  respect  to  cost)  this  change  in  the  Newsboy  problem 
should  occur.   Additionally,  we  wish  to  further  examine  the 
question  of  which  of  the  candidate  estimators  to  use. 

However,  it  was  hoped  that  the  estimators  of  X   or  S 
»  **  p     o 

would  come  close  to  having  desired  order  statistics  proper- 
ties, being  exact  for  the  case  of  the  uniform  distribution 
(equations  (9),  (10),  (11)  and  (12)).   For  example,  for  the 

r   order  statistic  X,  ,  it  has  been  shown  that 

(r) 

Fn(E[Xr    ^])  =  — ~r   for  a  uniform  distribution,  but 
D      u  (r J     n+1 

FD(E[X(r)])  <  -Jj.  for  convex,  and  FD(E[X(r)])  >  ^  for 

concave  cumulative  distribution  functions  of  demand.  [4]. 

Unfortunately,  there  is  not  an  analytical  measure  to  show 

how  much  the  slack  in  these  inequalities  would  change 

depending  upon  the  amount  of  convexity  or  concavity  of  the 

demand  distributions. 

Also,  for  the  third  candidate  estimator  of  X  ,  some 

P 

additional  assumptions  were  made  that  might  not  hold  for 
every  distribution. 
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In  summary,  our  reasons  for  simulation  are: 

1.  To  find  the  proper  time  (period)  to  begin  the 
transition  phase  (to  use  one  of  the  candidate  estimators) , 
leaving  the  "Uncertainty"  case  and  the  minimax  decision  rule, 

2.  To  provide  information  to  supplement  incomplete 
analytical  knowledge  of  order  statistics  regarding  popula- 
tion quantiles, 

3.  To  evaluate  the  impact  of  the  approximations  made 
in  developing  the  third  candidate  estimator,  X,  »,  over 
various  demand  distributions. 

B.   DESIGN  FOR  SIMULATION 

The  simulation  program  was  written  in  the  FORTRAN  IV 
language.   In  every  simulation  run,  inventory  stocking 
decisions  are  made  successively  using  one  particular  demand 

distribution  and  one  particular  value  of  the  quantile 

Cso 
(P  =  — p — t-j=s —  ) ,  say,  demands  generated  from  a  normal  parent 

c   so 

for  quantile  0.3  ...  etc.   The  demand  distributions  and 

various  quantile  values  used  in  the  simulation  are  given  in 

the  next  section.   Each  run  consisted  of  50  successive 

decision  periods  starting  with  no  demand  information  other 

than  the  range  estimate,  and  replicated  60  times. 

As  output  values  of  the  program,  costs  for  each  period 

were  averaged  over  the  number  of  replications.   These  average 

costs  per  period  were  calculated  for  five  different  cases: 

1.  The  Uncertainty  case,  with  the  minimax  decision  rule 
being  used,  denoted  by  M. 

2.  The  ideal  case  (Risk),  with  the  quantile  value  known, 
denoted  by  ID. 
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3.  The  transition  phase  case  using  estimator  X ,  ,,    denoted 
by  R.  Lrj 

4.  The_transition  phase  case  using  estimator  X,  ,  ,  denoted 

by  R.  l 

The  minimax  (for  Uncertainty)  and  ideal  (for  Risk)  cases 
would  give  constant  inventory  levels,  i.e.,  the  same  amounts 
all  the  way  through  the  simulation  run.   For  Uncertainty,  the 
range  of  demand  was  approximately  known  and  would  not  change 
unless  it  was  revised,  while  for  Risk  conditions  it  was 
assumed  that  the  distribution  of  demand,  and  thus  the  quan- 
tile,  was  known. 

Transition  phase  optimum  inventory  level  values  generated 

by  X,  -.  ,  X,  -.  and  X,  .  were  subject  to  change  from  period  to 
}       (r)    (r)       (r) 

period  until  finally  they  converge  to  the  value  associated 

with  the  Ideal  case.   For  the  first  two  periods,  X,  ^,  X,  ., 

v  (r)    (r) 

X-  ,  would  not  lead  to  an  inventory  level,  since  they  need 
at  least  two  data  points. 

The  program  was  run  for  50  periods,  and  replicated  60 
times.   So,  period  costs  for  each  transition  phase  procedure 
were  averaged  over  60  replications.   Program  steps  can  be 
ordered  in  the  following  way: 

1.  Define  the  demand  distribution  yielding  the  quantile 
p,  give  the  stockout  cost  C  ,  the  carrying  cost  C  ,  and  the 
range  of  demand. 

2.  Compute  and  set  minimax  and  ideal  case  optimum 
inventory  values. 

3.  Generate  a  random  demand  value  from  the  specified 
demand  distribution. 
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4.  Find  the  single-period  cost  for  each  case  (minimax, 

Xr  ->  ,  X,  >,  X,  ^  and  Ideal)  by  comparing  their  inventory 

levels  and  the  demand  generated  for  this  period.   Record 

these  costs.   Find  the  next  period  inventory  levels  revising 

X,-  ^  ,  Xr    -,  ,  X,   .  and  go  to  the  next  period, 
(r)'   (r)    (r)     5  F 

5.  Repeat  steps  3  and  4  for  50  periods. 

6.  Again,  beginning  from  the  first  period,  repeat 
step  5  for  60  times  (replications) ,  and  at  the  end  of  60 
replications  find  the  average  costs,  and  the  sample  variances 
of  costs  over  60  replications,  for  each  period. 

7.  Tabulate  the  results. 

The  results  obtained  from  this  simulation  will  be  examined 
in  part  D. 


C.   SIMULATION  PARAMETERS 

In  this  simulation  a  variety  of  possible  forms  of  distri- 
butions, including  continuous,  discrete,  symmetric,  skewed 
(long- tailed) ,  etc.,  were  chosen.   Demand  distributions  and 
their  parameters  used  in  this  simulation  were  the  following: 


Distribution 

Uniform 

Poisson 

Chi-square 

Lognormal 

Normal 


Parameters 

the  range  over  [0,30] 

with  mean  10  units 

with  10  degrees  of  freedom 

generated  by  a  normal  parent 
having  a  mean  of  3  and  a  variance 
of  1 

with  mean  35  and  variance  100 
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For  D    in  the  Minimax  procedure,  u  +  3t  is  used  for 
max  ^  r 

so 
tailed  distributions.   Quantiles  p  =  ^ — ^-= —  ,  tested  were 

c   so 
0.10,  0.30,  0.50,  0.70,  0.90,  0.95  for  each  distribution. 

These  were  determined  so  that  C   +  C  was  constant  (in  this 

simulation  C   +  C   =4.0  was  used).   So,  for  each  {quantile, 
so    c 

distribution}  combination  the  differences  of  the  expected 
period  costs  would  come  from  the  performance  of  the  estimators 
of  X  ,  as  was  shown  earlier  by  equation  (9) . 

In  the  next  section  we  shall  present  and  discuss  some 
of  the  results  from  the  simulation. 


D.   SIMULATION  RESULTS 

In  this  part,  simulation  results  will  be  evaluated  and, 
for  the  Newsboy  problem,  a  transition  phase  procedure  will 
be  proposed.   In  presenting  these  results  we  will  identify 
the  quantile  value  as  p.   Tables  I  and  II  give  overall  mean 
costs  and  standard  deviations  of  cost.   Then,  Table  III  will 
show  at  what  period  the  transition  procedure  should  begin 
for  a  particular  distribution,  and  quantile.   Table  IV  will 
give  relative  efficiency  of  the  best  transition  procedure 
to  the  Ideal  case  (decision  under  Risk) . 

On  the  next  page,  Table  I  gives  the  overall  mean  values 
of  the  cost  per  period  (over  50  periods  and  60  replications) 

for  each  transition  case,  using  C   +  C    4.00.   As  can  be 

'     °      so    c 

seen  from  Table  I,  there  does  not  seem  a  significant  differ- 
ence using  R,  R,  or  R,  even  though  R  looks  quite  consistent 
at  the  tail  (  p=  0.95).   One  reason  may  be  the  fact  that  the 
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TABLE  I 


Average  Costs  Over  50  Periods  and  60  Replications 
having  Cc  +  Cso  =  4.0 

M=Minimax,  P^X..,  R=X,r.  ,  I=X ,.  ,  ID=Ideal  Cases 


p 

Uniform 

Poisson 

Chi-Square 

Lognormal 

Normal 

M 

5.49 

3.23 

3.01 

67.00 

11.66 

R 

6.35 

2.38 

2.91 

14.50 

8.09 

0.1 

R 

6.31 

2.38 

2.88 

14.50 

8.09 

R 

6.30 

2.38 

2.87 

14.48 

8.07 

ID 

5.49 

2.04 

2.56 

11.44 

6.93 

M 

12.68 

4.60 

5.92 

255.96 

20.04 

R 

13.58 

4.55 

5.98 

36.34 

14.25 

0.3 

R 

13.55 

4.52 

5.96 

36.23 

14.18 

E 

13.46 

4.56 

5.94 

36.38 

14.16 

ID 

12.68 

4.26 

5.62 

30.15 

13.17 

M 

15.08 

5.97 

11.08 

339.86 

15.71 

R 

16.03 

5.35 

7.47 

51.22 

15.80 

0.5 

R 

15.97 

5.33 

7.45 

51.37 

15.78 

R 

15.91 

5.34 

7.45 

51.35 

15.73 

ID 

15.08 

4.98 

6.98 

43.49 

14.80 

M 

12.75 

7.33 

11.81 

298.67 

14.97 

R 

13.75 

4.74 

7.21 

57.24 

14.24 

0.7 

R 

13.66 

4.74 

7.21 

57.23 

14.19 

R" 

13.60 

4.72 

7.20 

57.66 

14.19 

ID 

12.75 

4.44 

6.67 

48.72 

13.24 

M 

5.39 

4.37 

6.26 

132.42 

9.33 

R 

6.31 

2.68 

4.40 

42.36 

7.76 

0.9 

R 

6.30 

2.67 

4.38 

42.38 

7.76 

R 

6.30 

2.66 

4.36 

42.58 

7.73 

ID 

5.39 

2.33 

3.71 

35.91 

6.70 

M 

2.83 

2.38 

3.34 

70.42 

5.27 

R 

3.77 

1.75 

3.00 

31.80 

5.04 

0.95 

R 

3.76 

1.75 

2.99 

31.74 

5.02 

R" 

3.76 

1.75 

3.00 

32.44 

5.04 

ID 

2.83 

1.41 

2.27 

25.47 

3.91 
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difference  between  variance  of  X,  -.  and  variance  of  X  ,  -.  , 
was  not  much  (Figure  3) .   But  if  one  wishes  to  see  the  dif- 
ference no  matter  how  much  it  is,  the  following  normalization 
by  ideal  costs  may  be  considered: 

1.  For  a  distribution  from  Table  I,  for  every  quantile 
R.-ID 

p_  find  -y=j —  ,  and  sum  this  value  over  all  quantiles  g_  for 

every  R. ,  where  R-  is  R,  R  or  R. 

2.  Repeat  step  1  for  every  distribution. 

3.  Pick  the  R.'s  which  have  the  smallest  sum  among 
R,  R,  R  for  every  distribution  as  their  best  procedure. 

The  results  of  this  normalization  were  as  follows: 

Uniform  Poisson  Chi-square  Lognormal  Normal 

R    0.871  0.767  0.859  1.253  0.839 

R    0.845  0.752  0.831  1.251  0.824 

I   0.828  0.755  0.821  1.295  0.816 

The  best  transition  procedures  were  to  use  R  for  uniform, 
R  for  poisson,  R  for  chi-square,  R  for  lognormal,  and  R  for 
normal  distributions.   To  see  the  effects  of  the  very  right 
tail  of  a  distribution  on  transition  procedures,  p  =  0.99 
was  used  for  normal  and  chi-square  distributions.   At  this 
quantile  the  minimax  costs  were  below  R,  R,  and  R  costs. 
This  says  minimax  can  be  used  at  the  very  tail  of  a  distri- 
bution along  the  transition  phase,  providing  very  accurate 
prediction  of  the  range  of  demand  (say  +51).   Otherwise,  if 
we  remember  from  Table  I  for  p  =  0.95,  R,  R  and  R  were 
better  than  minimax  approach. 
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From  Table  II  (overall  standard  deviations  over  50  periods 
and  60  replications) ,  it  is  seen  again  that  there  is  not  much 
difference  of  standard  deviations,  even  though  R  has  a  better 
trend  (smaller  sample  variance  of  cost)  among  .R,  R  and  R. 

From  these  two  tables  it  can  be  said  that  there  is  not 
much  difference  of  the  costs  yielded  using  R,  R  and  R,  but 
if  an  order  is  necessary  it  should  be  the  one  using  R  for  p 
less  than  0.95,  and  R  at  the  tail  (i.e.,  when  p  is  greater 
than  0.95) . 

During  the  simulation,  as  the  number  of  periods  increase 

the  estimators  converge  to  X  or  the  true  S  ,  while  the 

minimax  and  ideal  case  inventory  levels  were  constant. 

Since  the  values  of  X  ,-  -,  ,  X,  •,  ,  X,-  .  usually  change  from 

(r)    (r)    (r) 

period  to  period  because  they  are  updated  by  the  most  recent 
demand  value,  they  will  generally  be  associated  with  a  larger 
variance  of  cost  than  the  minimix  and  ideal  cases.   This  is 
reflected  in  Table  II.   Sometimes  the  minimax  procedure  had 
a  smaller  sample  variance  of  cost  (while  possibly  yielding 
the  largest  average  cost)  than  the  ideal  case.   This  may  be 
explained  as  follows: 

The  ideal  solution  we  are  using  relates  to  minimizing 
expected  cost,  not  the  variance  of  cost.   These  two  minimi- 
zation points  are  generally  different,  depending  on  the 
distribution  function.   If  the  minimax  solution  is  different 
from  the  ideal  (expected  cost)  inventory  level  and  closer  to 
the  minimum  variance  solution,  it  would  yield  a  smaller 
variance  of  cost.   Further  derivations  about  the  variance 
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TABLE  II 

Sample  Standard  Deviations  over  50  Periods 
and  60  Replications 

M=Minimax,  R=Xfr-.,  R=Xfr.  ,  R^X,.  ,  ID=Ideal  Cases 


p 

Uniform 

Poisson 

Chi-Square 

Lognormal 

Normal 

M 

3.15 

1.26 

1.79 

43.46 

3.77 

R 

6.36 

2.69 

3.27 

27.46 

9.89 

0.1 

R 

6.26 

2.65 

3.21 

27.47 

9.70 

R 

6.26 

2.59 

3.18 

27.47 

9.53 

ID 

3.15 

1.62 

2.26 

15.90 

6.62 

M 

7.46 

3.18 

4.61 

76.10 

10.33 

R 

9.00 

3.60 

4.93 

57.64 

12.17 

0.3 

R 

8.84 

3.53 

4.89 

57.69 

11.97 

R 

8.85 

3.51 

4.86 

57.71 

11.99 

ID 

7.46 

3.16 

4.42 

45.93 

10.95 

M 

8.82 

4.30 

6.35 

66.95 

12.34 

R 

10.31 

4.14 

6.37 

84.29 

12.73 

0.  5 

R 

10.  23 

4.12 

6.33 

83.91 

12.68 

R 

10.23 

4.14 

6.44 

84.23 

12.75 

ID 

8.82 

3.84 

5.80 

72.72 

12.09 

M 

7.30 

3.47 

4.87 

46.79 

9.  62 

R 

9.08 

3.95 

6.66 

99.68 

11.65 

0.7 

R 

9.01 

4.04 

6.83 

99.68 

11.65 

R 

8.87 

4.09 

6.76 

99.23 

11.67 

ID 

7.30 

3.55 

5.66 

92.87 

10.42 

M 

3.16 

1.25 

1.87 

32.45 

3.90 

R 

6.59 

3.28 

5.72 

99.99 

9.10 

0.9 

R 

6.55 

3.33 

5.76 

99.23 

9.09 

R 

6.55 

3.29 

5.69 

98.05 

9.01 

ID 

3.19 

2.57 

4.10 

94.42 

6.34 

M 

1.66 

0.69 

1.12 

31.66 

2.07 

R 

6.50 

2.91 

5.36 

94.64 

8.50 

0.95 

R 

6.50 

2.95 

5.34 

92.85 

8.47 

R 

6.50 

2.94 

5.23 

92.74 

8.46 

ID 

1.66 

1.49 

3.15 

84.97 

4.11 
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minimization  solutions  of  the  Newsboy  problem  are  given  in 

Appendix  B. 

The  simulation  results  about  the  period  to  begin  a 

transition  phase  are  given  in  Table  III.   Shown  are  the 

first  periods  in  which  the  average  period  cost  (over  60 

replications)  using  the  estimator  was  below  that  obtained 

using  the  minimax  rule.   In  the  simulation  it  was  observed 

that  once  costs  using  an  estimator  fell  below  the  minimax 

costs,  they  stayed  below  for  the  remaining  time  periods. 

Simulation  results  showed  that,  on  the  average,  at  the  fifth 

period  we  were  supposed  to  begin  a  transition  procedure, 

ceasing  to  use  the  minimax  procedure.   In  other  words,  for 

the  first  five  periods  the  large  variability  of  candidate 

estimators  of  S   (or  X  )  caused  a  larger  average  cost  than 
op 

the  minimax  procedure  yielded.   As  a  very  special  case,  if 
the  distribution  of  demand  is  uniform,  minimax  and  Ideal 
case  solutions  would  match  each  other  and  thus  the  uniform 

distribution  is  omitted  from  Table  III.   Also  by  coincidence, 

so 
if  p  =  p — ttt —  is  at  or  close  to  an  intersection  of  the 
Lc  Lso 

implied  uniform  CDF  of  minimax  procedure  and  the  actual 
distribution  function  of  demand,  then  the  minimax  procedure 
may  yield  the  Ideal  case  solution.   Otherwise,  minimax  in- 
ventory level  has  always  a  higher  cost  than  the  Ideal  case. 
In  Table  III  the  large  numbers  are  those  specific  values 
which  relate  to  a  quantile  P  that  is  close  to  minimax-- 
actual  CDF  intersection  point.   So  it  takes  time  for  R,  R 
and  R  to  yield  costs  less  than  the  minimax.   No  matter  how 
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TABLE  III 

Periods  to  Begin  Transition  Phase 
Leaving  the  Minimax  Procedure 


Quantile 

(P) 

Poisson 

Chi-Square 

Log: 

normal 

Normal 

0.1 

6 

6 

5 

6 

0.3 

10* 

10* 

4 

4 

0.5 

6 

4 

4 

10* 

0.7 

4 

4 

4 

6 

0.9 

5 

5 

4 

5 

0.95 

6 

5 

4 

6 

*  Ideal-Minimax  intercept  points  (or,  close  to  these  points) 
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long  it  will  take  to  fall  below  the  minimax  cost  (here, 
nearly  at  period  ten)  from  simulation  outputs  it  was  observed 
that  after  the  fifth  or  sixth  period  the  cost  difference 
between  the  minimax  and  the  transition  procedures  was  not 
large.   This  also  can  be  predicted  from  the  pattern  values 
of  Table  III.   The  numbers  in  Table  III  are  relevant  for  all 
candidate  estimators,  i.e.,  R,  R  and  R  all  fell  below  the 
minimax  cost  simultaneously  at  the  same  period. 

As  a  measure  of  efficiency  for  a  transition  procedure, 
the  ratio  of  average  cost  induced  by  this  procedure  to  Ideal 
average  cost  is  used.   Table  IV  gives  these  ratios  for  the 
average  costs  over  second  and  fifth  10  periods  for  R  (since 
R  seems  the  promising  transition  procedure  among  the  other 
ones)  . 

We  will  examine  the  data  of  Table  IV  for  the  second  10 
period,  assuming  that  the  minimax  has  been  left  and  the 
transition  phase  has  begun.   We  also  believe  that  convergence 

toward  the  true  value  of  Ideal  inventory  level  S   for  the 

J  o 

third,  fourth,  ...  10  periods  will  not  be  worse  than  the 
convergence  accomplished  during  the  second  10  period. 

If  the  transition .procedure  (here,  R)  was  operated  at 
the  tail  (i.e.,  p  >_  0.9),  then  the  efficiency  ratio  for  the 
second  10  periods  was  on  the  average  1.16  and  more  (worse). 
As  the  quantile  gets  close  to  1.0,  the  minimax  procedure  was 
yielding  much  better  results  than  a  transition  procedure 
(assuming  the  range  of  demand  is  predicted  quite  well) . 
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TABLE  IV 

The  Ratio  of  the  Second  and  Fifth  Ten-Period 
(11  to  20,  41  to  40)   Average  Costs   lo  the 
Ideal  Case  (Under  Risk)  Costs,  Using  X ,  ^ 
Transition  Procedure        ^  ' 

Quantile  (p)   Uniform  Poisson  Chi-Square   Lognormal   Normal 

S     1.063     1.065      1.058       1.033     1.108 
0.1 

F     1.000     1.023      1.046       1.002     1.008 


0.3 

S 
F 

1.065 
1.007 

1.059 
1.041 

1.050 
1.003 

1.040 
1.004 

1.066 
1.008 

0.5 

S 

F 

1.056 
1.013 

1.069 
1.021 

1.039 
1.013 

1.041 
1.004 

1.077 
1.007 

0.7 

S 

F 

1.073 
1.000 

1.069 
1.001 

1.061 
1.006 

1.057 
1.006 

1.054 
1.004 

0.9 

S 
F 

1.102 
1.008 

1.113 
1.016 

1.128 
1.052 

1.190 
1.026 

1.096 
1.027 

0.95 

S 
F 

1.170 
1.002 

1.193 
1.040 

1.211 
1.054 

1.240 
1.062 

1.183 
1.055 

Notes:  . 

S  stands  for  "SECOND,"  F  stands  for  "FIFTH." 

When  p  is  less  than  0.9  overall  ratio  is  1.060  for  second 
and  1.010  for  fifth  10  periods. 
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For  the  rest  of  the  quantiles  where  p  is  less  than  0.9, 
this  ratio  for  the  transition  procedure  was  almost  constant 
1.060.  To  test  the  effects  of  quantiles  and  the  effects  of 
the  form  of  demand  distributions  on  the  efficiency  ratio 
(1.060),  a  two-way  Analysis  of  Variance  could  be  applied  to 
the  data  of  Table  IV  for  quantiles  less  than  0.9. 

Different  distributions  were  significant  at  a  level  0.03, 

and  the  different  quantiles  were  significant  at  a  level  of 

0.01.   This  result  was  explained  that  for  C   +  C   =4.0 

r  so    c 

cost  units,  whatever  the  distribution  of  demand  was  and  for 
any  quantile  less  than  0.9,  the  ratio  of  the  cost  induced  by 
transition  phase  procedure  (X,  •>  ,  here)  to  the  Ideal  cost 
(expected  cost  under  Risk)  would  be  practically  constant  at 
1.060  for  second  10  period,  and  would  not  exceed  this  amount 
for  the  further  periods  (for  example,  at  the  fifth  10  period 
this  ratio  was  1.01). 

On  the  next  part,  as  the  summary  of  this  chapter,  and 
using  the  results  obtained  in  it,  a  Single-Period  inventory 
(Newsboy)  problem  transition  phase  procedure  is  suggested. 

E.   PROPOSED  NEWSBOY  TRANSITION  PROCEDURE  ALGORITHM 

(1)  Approximate  the  bounds  of  Demand  (Range) . 

(2)  Set  p  =  r   r —  . 
Lc  +  Lso 

(3)  If  p  >  0.90,  go  to  step  6,  otherwise  go  to  step  4. 

(4)  Use  minimax  for  the  first  5  or  6  periods. 

(5)  For  6  <   n  £  50,  use  R  (unless  a  reason  for  R,  R) , 
and  go  to  step  10. 

(6)  If  p  <  0.95,  go  to  step  8,  otherwise  go  to  step  7. 
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(7)  If  step  1  is  accurage,  may  use  minimax  for  all  periods 
and  go  to  step  10,  but  if  it  is  not  accurate,  go  to 
step  8. 

(8)  Use  minimax  for  the  first  5  or  6  periods,  then  go  to 
step  9. 

(9)  For  6  <_  n  <_   50,  use  R  and  go  to  10. 

(10)   Analyze  data. 

*May  fit  a  probability  distribution  to  the  data  and  begin 
Decision  Under  Risk  case. 


*May  use  R,  R  or  R  for  n  >  50  too.   Since  they  converge 
X  quite  fast,  the  risk  of  using  them  after  50  periods  hope- 
fully will  be  nearly  equal  to  the  risk  of  Ideal  case. 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 
FOR  FURTHER  STUDY 

For  the  Newsboy  problem,  a  nonparametric  decision  pro- 
cedure that  allows  transition  from  Uncertainty  to  decision 
under  Risk  has  been  developed.   Three  estimators  of  the 
demand  corresponding  to  the  p   quantile  were  obtained, 

where  p  =  p — —^ —  .   Testing  these  estimators  over  several 
Lc  +  Lso 

distributions  of  demand  with  Newsboy  cost  criteria,  it  was 
found  that  excluding  cases  such  as  drastic  skewness,  pooling 
the  (r-l)st  and  (r+l)st  order  statistics  for  R  provided 
essentially  the  best  estimator,  where  R  and  R  also  worked 
well. 

It  was  also  found  that  for  the  distributions  tested  and 

SO 

excluding  (•*-= —  r  )  values  very  close  to  0.0  or  1.0,  that 
Lso   Lc 

the  switching  Newsboy  decision  making  from  minimax  solutions 
to  quantile  estimator  could  be  made,  possibly  at  the  sixth 
period. 

For  further  work,  some  parametric  studies  such  as  log- 
normal  distribution,  or  beta  and  gamma  families,  depending 
upon  whether  the  population  of  demand  is  finite  or  not,  may 
be  tried.   Then,  given  that  the  decision  maker  does  not  know 
what  the  distribution  of  demand  is,  which  procedure,  nonpara- 
metric or  parametric,  gives  better  response  to  Uncertainty 
may  be  compared.   As  another  recommendation,  some  work  can 
be  afforded  to  improve  nonparametric  approaches  further. 
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APPENDIX  A 

ESTIMATION  FOR  POPULATION  QUANTILES 
BY  ORDER  STATISTICS 


In  this  Appendix  order  statistics,  especially  those 
concerning  population  quantiles,  will  be  briefly  reviewed. 
For  a  further  study  of  order  statistics,  references  [1], 
[2]  and  [3]  are  cited. 

The  Joint  Distribution  of  n  Order  Statistics 

Order  statistics  theory  is  mainly  based  on  continuous 
distribution  functions.    In  reference  [2]  for  discrete 
parents  the  joint  distribution  of  order  statistics  is  derived 
as  well. 

Let  X, ,X2,...,X  be  a  random  sample  from  a  probability 

distribution  function  Fv.   If  X, . . , i=l, 2 , . . . ,n  is  the  i 

X       (lj     '  '    ' 

largest  value  in  this  sample,  then  X,-,-.  <  X.2.  <  ...  <  X,  -. 
are  the  order  statistics  of  a  random  sample  X, ,X2,...,X  . 
Here,  X,  ■.    is  called  the  r   order  statistic.   Since  FY  is 
a  distribution  function,  by  the  probability  integral  trans- 
formation U  =  F„(x)  will  have  uniform  distribution  over  [0,1], 
i.e.  , 

!u,  for  0  <  u  <  1 
0,  otherwise. 

Letting  F„(x)  =  u  with  F„  (X)  <  u  if  and  only  if  X  <_  x, 
the  following  equations  hold:  ,-,,. 

FyCu)  =  P(U<u)  =  P{Fx(X)<u}  =  P[X<Fx_1(u)]  =  Fx[Fx1(u)]  =u 
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Since  FY(X)  is  a  monotone  increasing  function  and  by 


equation  (16)  there  is  a  one-to-one  correspondence  between 
FX(XC1))  <  FX(XC2))  <  ...<Fx(X(n))  andX(1)  <  X(2)  ...  < 
Let  fy   be  the  density  of  FY,  so  that 


(n) 


fY   Y      Y  (x,,x7,...,x  )  is  the  joint  density  of 
1  L  n 

the  sample  X,  ,X2»  •  •  •  »X  ..   If  each  X.  ,   i=l ,  2  , .  .  .  ,n 

is  independent  and  identically  distributed,  then 

n 
f, 


X 


Y  Y   (^  "|  J-^?  >  •  •  •  >^n  J  ^tY^X-J 

i >A2>  *  •  '  '  n  i=l 


Since  there  exists  n!  permutation  of  those  sample  values, 

n 

fY    v        y    (.Y;,y2>---  >YJ    =  n!  =  f  (y  )  ,    (17) 
fnf21'""*fn*) 

is  the  joint  density  of  n  order  statistics. 

From  the  previous  references  cited: 

n!  r-1         n-r 

£X(r)Cyr3  =  (r-1)!  (n-r)!  fFX^r)]   ^"V^    W' 

(18) 

When  F„  is  a  uniform  distribution  over  (0,1),  this  equation 

turns  to  be, 

£X(r)^  =  (r-1)!  (n-r)!  ™       ^   '  >  ™ 

and  the  following  sequence  of  equations  are  given, 
E[Xr   *]  =  ^.?,',fTt.T.s,  /  u1"^"1  (l-uf"rdu 


(r)  J    (r-1)!  (n-r)!  f) 


and 
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k     (r+k-1) (r+k-2) (r+l)r 

E[X(r)  ]  =  (n+k)(n+k-l) 0+2)  (n+1)  '  (20) 

where  k  is  a  positive  integer.   For  k=l: 

E[x(r)]  -  &  .  en) 

By  equations  (20)  and  (21) 

Var[XCr)]    -  E[X(r)2]    -    [E(X(r))]2  -    (^^)2     .      (22) 

and, 

E[Xf         -    xr    J    =   -il±^ ,  (23) 

L    Cr)  (s)J         (n+l)(n+2) 


From   equations    (6)    and    (8), 

Cov[X(r),X(s)]    =    E[X(r)X(s)]    -    E[X(r)]E[X(s)],    and 

Cov    [XX]    =        r(n-s+l)         •  (24) 

(r)       (sj  (n+l)Z  (n+2) 


Now,  we  wish  a  transition  from  the  uniform  distribution 
to  a  general  distribution  F„.   Using  the  relations  in 
equation  (16) , 

U(r)  =  Fx(X(r))  ,     or    X(r)  =  F,"1  (U(r))  , 

where  U,  ^  is  the  r   order  statistic  from  Uniform  (0,1). 

Then,  E[X(r)]  =  E^"1  (U(r))]  *  F^1  [E(U(r))j  =  F^1  (^) 

(25) 
Furthermore,  the' references  given  at  the  beginning  of  this 
appendix  show  that,  for  continuous  distributions  having 
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n.  =  [np.  +  l] ,  and  0  <  p1 <  p2   . . .  <  p^  <  1  , 

the  joint  distribution  of  X,   i>Xfn  V'*'X(nJU  asy7ntotically 
tends  to  an  ^-dimensional  normal  distribution  with  means 
p..  ,  p-,...,  p*  ,  and  covariances, 


Cov[X(n.)'X(nv)]  =  nf  (X  )  f  (X  )  '  Provided  i   1 

v  l7    v  KJ  Pi       Pic 


k 


(26) 


So  the  variance  of  X,   .  will  be, 


PiCi-Pi) 
Var[Xr„  0  = 


Cni5     n{f(Xp  )}2 
Having  n.=  r,  it  can  be  proven  that, 

^-V     ^  (XCr)'XPi} 
U  =     1        =  -  fx(Xn  )  %  N(0,1)    (27) 

/Var(X(n  })    /p^l-p.)      *  Pi 

Equation  (27)  is  quite  important  for  this  study.   It 

says  that,  whatever  the  distribution  function  is,  the  r 

order  statistic,  X^  -.  ,  is  approximately  normally  distributed 

around  X   (distribution  free  property) .   This  also  means 

that  asymtotically  Xr    >,    is  an  unbiased  estimator  of  X  . 

3       (r)  p 

Let  Gn(X)  be  the  empirical  distribution  function  of 
order  statistics  X., -.  ,X  ,2^  , .  .  .  ,X,  >   having, 

0  ,  if  X  <  X,^ 
GD(X)  =    \       j  ,  if   X(j)  <  X  <  XCj  +  1) 

1  ,  if  X  >  x(n) 
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In  Reference  [3]  it  is  said,  by  the  Glivenko-Centelli 
theorem  that,  Gn(X)  converges  uniformly  to  F„(X)  and  nGD(X) 
is  binomial  (n,Fx(X)),  i.e., 

p[nGD(X)  =  j]  =  p[GD(X)  -  1]  =  (")[Fx(X)]J[l-Fx(X)]n-J, 

(.28) 

Point  and  Interval  Estimation  of  Population  Quantiles 

For  point  estimation  of  a  population  quantile,  we  assume 

that  X   is  the  p   quantile  of  the  distribution  function  Fv, 
p  ^  X 

and  P(X<Xp)  =  FxCXp)  =  p  ,  or   F^1  (p)  =  Xp . 

Then,  by  equation  (27)  the  r   order  statistic  is  a  consistent 
and  unbiased  estimator  of  the  p   quantile,  where 


r   np     ,   if  np  is  integer 
r  ■  •  C 

„  [np+1]  ,   otherwise  ,   and  p  =  ^ — ^ 


so 

+C 
so   c 


For  interval  estimation  of  population  quantiles,  we  are 
looking  for  r  and  s  such  that, 

PfX,  ■>  <  X   <  X,  n)  ■  y  ,  where  y  is  percent  confidence, 
(r)    P     (sj 

Now, 

PtX(r)  <  Xp  <  X(s)^  =  PCX(r)  <  V  '  P(X(s)<  V'  °r 
by  equation  (16) 

P(X(r)  <   Xp  <  X(s)^  =  PCFX<X(r))  <  P)  "  P(VX(s))  <  P3- 

Since  p  =  FY(X  ),  by  equation  (28) 
a  p 

P(xf  <x  <x   )  -  I   C^P^i-pf"1-  s  C?)  p^i-pf"1 

i=r  i=s 

s-1  „ 
-   Z  Ci)p1d-pf"1-  Y  (29) 
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To  summarize  these  derivations: 

1.  For  large  n  the  r  order  statistic  is  a  consistent, 
unbiased  point  estimator  of  the  p  quantile  for  any  distri- 
bution. 

"f"  li 

2.  Confidence  intervals  for  the  p  quantile  as  given 
by  equation  (29) ,  distributed  by  the  Binomial  distribution 
(n,p). 

3.  As  n  gets  larger,  X,  .  converges  to  X   and  also  the 
empirical  distribution  function  Gp(X)  converges  to  actual 
cumulative  distribution  function  Fv(X) . 
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APPENDIX  B 

EXPECTED  COST  AND  VARIANCE  MINIMIZATION  POINTS 
OF  THE  NEWSBOY  PROBLEM 


In  this  Appendix  we  will  attempt  to  explain  the  differ- 
ence between  minimization  of  the  expected  cost,  and  mini- 
mization of  the  variance  of  the  cost,  for  the  Newsboy 
problem. 

Decision  makers'  utility  functions  are  involved  with 
both  expected  cost  (or  return) ,  and  the  variance  of  this 
cost.   If  a  decision  maker  is  risk  prone,  the  variance  of 
the  cost  may  not  affect  him,  and  he  just  tries  to  minimize 
expected  cost.   On  the  other  hand,  if  he  is  risk  adverse, 
then  he  may  give  up  some  of  the  cost  minimization  value 
for  a  smaller  variance  (a  trade-off  between  expected  cost 
and  variance  of  cost).   As  an  example  [5],  let 

2 
U(x)  =  x  -  bx  ,  be  a  utility  function.   Then, 

E[U(x)]  =  E[x]  -bE[x2]  =  E[x]  -  b(Var(x)  +  [E(x)]2). 

If  we  try  to  optimize  this  value,  we  not  only  deal  with 
expected  value  of  x,  but  also  variance  of  x. 

On  the  following  pages,  the  variance  minimization  points 
of  some  distributions  related  with  the  Newsboy  problem  are 
derived. 
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Let  VarCC)    =    E[C2]    -   E[C]2    ,    and 


E[C]    =   Cc/    CS-x)£D(x)dx   +    Cso/(x-S)fD(n)dx 
0  j 

s  °° 

E[C2]    =   Cc2/    CS-x)2£D(n)dx   +   C^/ (x-S)  2f  D(n)  dx  . 

U  o 


To  minimize  Var(C),  with  respect  to  S: 


dVar(C)         n      . 

— ,  ^—■'    =    0,    is   necessary. 


So, 


dVar(C)    =   dE[CZ]    _    ,    f    ,    dE[Cj    m   Q 
d~S  dT~        ztlLJ       dS  u' 


(30) 


Then, 


2  2  2   °° 

^LCJ.  .    2C2f    (S-x)£D(x)    -    2CsQ   /Cx-S)£D(x)dx 

U  o 


2E 


[c]    dEXCl=    .    2£[c]     {  £c(x)dx-Cso7£D(x)dx} 


=    -    2E[C]{CcFD(S)  -  Cso(l-FD(S))} 


dE[CT] 


dE[C] 


Since,   ^y-  =    2E[C]    ^ 


s=s. 


s=s. 


then   it   yields 


C1(SF(S)  -  /xfn(x)dx}    -    cM/xfn(x)dx  -  S(l-Fn(S))}    = 


0 


■D 


so    c        D 


D' 


(CC{SF(S)  -    /xfD(x)dx}    +    Cso{/x£D(x)dx  -  S(1-FD(S))}) 

U  o 


(Fn(S)  (C     +  C      )  -  C      ) 
^Dc        so  so 


(31) 
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When  this  is  applied  to  a  uniform  distribution  over 

CO,  D   ),  i.e.,  fn(x)  =  — !—  ,  0  <  x  <  D     and  Fn(S)  =  =S 
*■  ■  '   maxJ  '      '   D^  J        D '    —   —  max        D*-  J        D 


max  "max 


then,  the  left  hand  side  (LHS)  of  equation  gives, 

LHS  =  S2(C  2  -  C  2)  ^ —  +  C  2  S  -  C  2  -5S£  , 

v  c     so^  2D  so      so   2 
-  max 


Similarly,  right  hand  side  (RHS)  of  it  yields, 

C  +C  q2  CD 

RHS  =  (  Ctt — — )  S  -  C   )  (J? (C  +  C   )  -  C   S  +  s°  max)  . 

v VD     ^     so; V2D     v  c   so;    so     2    ^ 
max  max 

After  we  equate  LHS  -  RHS  and  opened  them,  it  was  found  that, 

C   D 
S   =  so  max ,  (32) 

0    C   +  C 

c    so 

is  the  point  that  minimizes  the  variance  of  the  cost  for 
the  uniform  distribution.   This  point  is  exactly  the  same 
point  that  minimizes  expected  cost. 

Let  us  apply  the  same  procedure  to  a  bell-shape,  say, 
Beta  (a  =  2,  3  =  2) .   Then 

£D(X)  =  rja)f(e)  ^^  Cl'X)3"1  -  6x(l-x)  for  0  <  x  <  1, 

Fn(x)  =   /  6u(l-u)du  =  3x2  -  2x3,  and 
u  0 

1  -  Fjj(x)  ■  1  -  3x2  +  2x3,   /"  u  fD(u)du  =  ^(l-4x3  +  3x4)  , 

A. 


/uf.  (u)du  =  i  (4x3  -  3x4  ). 
0 

Putting  those  values  in  equation  (2),  gives  a  5  degree 

polynomial : 
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<Cc  +  Cso^2  s5  -  I(Cc  +  Cso^2  S'  +  3  <Cc  +  Cso)2  s3  + 


+  ifC  2  +  C  2  +  5C  C   )S2  -  (4C  2  +  5C  C   )S  +  ^(C  C   +  2C  2  )  =  0 
2^0   so     c  scr     v   so    c  so;    2  *•  c  so    so; 

(33) 
It  may  be  difficult  to  solve  this  equation,  but  it  is 

not  difficult  to  see  that  it  doesn't  have  the  same  solution 

C 
with   so    _  ,c  2   ,c  3       ,c  2   ,,c  3     _  n 

n    .r =  3S„  -  2S„  ,  or   3S„  -  2S„  -  p  =  0, 

C  +C       o     o  '       o     o   r     ' 
c  so 

which  is  the  minimum  expected  cost  solution. 
In  the  exponential  {A}  case, 

fD(x)  =  Ae~Xx,  FD(x)  =  l-e~Xx,  l-FD(x)  =  e~Xx   and 

f   u  f(u)du  =  f[l-e~Xx(l+  x)], 
0  A 

/°°u  f(u)du  =  ^[e"Xx(l+Ax.)]. 
x  A 

Again,  inserting  these  values  into  equation  (2)  gives, 


FD(S)    CcX 

s    U^c   '  ^34^ 


where  FD(S)  =  1-e 


c   so 

-AS 

CI 


If  we  compare  C34)  with  FD(SQ)  =  ^-|g =  l-e~XSo 

r  c   so 

1       c 
or  S   =  -  yln  (r   +f — )  >  we  will  see  that  they  are  different 

c   so 
points.   If  we  examine  equation  (34)  further,  it  can  be 

shown  that  S  >_   S  ,  i.e.,  the  variance  minimization  point  S 

is  closer  to  the  tail  than  the  minimum  expected  cost  point  S 
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